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We explore analytically the structure of relativistic shock and solitary wave solutions in collision- 
less plasmas. In the wave frame of reference, a cold plasma is flowing from one end and impacting on 
a low velocity plasma. First we show that under astrophysical conditions, a cold electron-positron 
plasma is unstable with respect to a two-stream instability in the interface between these regions. 
The instability heats the inflowing cold plasma rapidly, on a timescale comparable to the inverse of 
its plasma frequency. We then derive time-independent equations to describe the resulting hot state 
of the pair plasma, and describe the conditions under which the spatially uniform solution is the 
unique stable solution for the post shock conditions. We also examine plasmas composed of cold 
protons and hot electrons, and show that the spatially uniform solution is the unique stable solution 
there as well. We state the shock jump conditions which connect a cold, electron-proton plasma to 
a hot electron-proton plasma. The generic feature evident in all of these models is that the plasma's 
initial, directed kinetic energy gets almost completely converted into heat. The magnetic field plays 
, the role of catalyst which can induce the plasma instability, but our solutions indicate that the 

macroscopic field only gets amplified by a factor of ~ 3 in the frame of the shock. 
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I. INTRODUCTION 



> 

Collisionless shock waves are common in astrophysical environments. They are a generic product of cosmic ex- 
plosions, such as supernovae jl],|j (which result in non-relativistic shocks), and gamma-ray bursts ||,^| (which yield 
highly relativistic shocks). Understanding their structure and dynamics is crucial for modeling pulsar magnetospheres 
H , active galactic nuclei Q , structure formation in the intergalactic medium , and a local phenomenon such as the 
Earth's bow shock in the solar wind JsJ . Despite the wide range of applications related to collisionless shocks, their 
theoretical modeling is still at a primitive state, and so astrophysicists often resort to collisional fluid equations in 
modeling extra-solar systems with no proper justification. 

Because of the disparate physics used to describe collisional and collisionless plasmas, it is important to have an 
independent model of collisionless shocks. In a collisional plasma, the inter-collision mean free path sets the scale for 
the shock transition layer. Frequent collisions generate a Maxwellian distribution of post-shock particle energies, and 
all particle species eventually come to share a common temperature [p|pO|]. However, shocks also occur in plasmas for 
which binary collisions are completely negligible 1 11 l^,[l3| . In this case, plasma instabilities act through macroscopic 



electromagnetic fields to bring about the shock wave, and possibly establish thermal equilibrium 

Several aspects of collisionless shocks have been explored in the literature. Steady, soliton solutions have been found 
mm for cold, magnetically dominated pair plasmas. But when the plasma kinetic energy is comparable to the Poynting 



energy, particles are magnetically reflected and execute looping motions. Simulations have suggested |jnj,|rj,|15[ that 
these looping motions are unstable, and that a shock front forms in the region of magnetic reflection. Downstream 
of the shock, the plasma kinetic energy was seen to be mainly converted into thermal energy. Since the simulations 
allowed particle motion in only two dimensions, it is unclear whether or not the resulting momentum distribution 
is isotropic. These simulations have been complemented by analytic studies |T6|]ll] , |l8| of the two-stream instability, 
which is operative whenever the fluid is composed of counter-moving streams. Since the aforementioned looping 
orbits present such a case, these analytical arguments suggest that the two-stream instability would dissipate the 
particles' ordered motion on time scales comparable to the inverse plasma frequency. As the instability is saturated, 
the particles' momentum distribution isotropizes, and the magnetic field receives significant amplification. 

In several astrophysical situations, and in particular the generation of gamma-ray burst afterglows, it is important 
to understand what happens when a relativistic shock wave impinges on a cold plasma carrying a small amount of 
energy in its magnetic field, compared to the kinetic energy of the shock. In §11 we will investigate the stability of 
a multi-stream configuration that typically results at such an interface. We will derive a simple dispersion relation, 
valid when the plasma frequency is much greater than the cyclotron frequency, which indicates that this configura- 
tion is unstable. This instability acts on a time scale comparable to that observed in simulations, and provides a 
promising driver for heating the gas and creating the necessary jump conditions of collisionless, relativistic shocks. 
The instability eventually saturates and converts the ordered bulk velocity of the particles into a three-dimensional 
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velocity distribution (thus, playing a role similar to collisions in a normal fluid). In §111, we will develop a stationary 
fluid model to describe the the post-shock plasma. After developing the appropriate equations, we will examine two 
cases in particular. First, we will investigate hot pair plasmas (§IV), and show that no stationary, continuous shock 
solutions exist. We will further derive the conditions under which oscillatory solutions might exist, and show that they 
are similar to those obtained earlier for cold plasmas In §V, we will investigate the case of an electron-proton 
plasma. Here, the two species have very different dynamical length scales, and so we will construct a simple model in 
which the electrons are hot, but the protons are cold. As a first approximation, we will average over the behaviour 
of the electrons, and calculate the resulting fields and the motion of the protons. Finally, we will consider a plasma 
in which both the electrons and the protons are hot. We will derived the jump conditions, and briefly discuss the 
possibilities for spatial structure. §VI will summarize the main conclusions of this work. 



II. THE INSTABILITY 



Previous studies have revealed that a generic feature of an initially cold plasma carrying relatively little energy in 



its magnetic field is the development of multiple streams |11 12] , |l3| , |l4| , |l8[ . The condition that there be little energy in 
the magnetic field can be stated in terms of the dimensionless ratio a = B 2 /8irj 2 nmc 2 <§; 1, where B is the plasma's 
self-magnetic field, 7 is the Lorentz factor corresponding to the mean fluid velocity, n is the proper number density, m 
is the mass per particle, and c is the speed of light. For simplicity, we suppose that the direction of the magnetic field 
is initially perpendicular to the mean fluid velocity. Then the definition of a indicates that it is invariant with respect 
to Lorentz boosts in the direction of the bulk flow. This fact allows for a simple calculation of a in the plasma's rest 
frame; for electrons in the interstellar medium, a ~ 5 x 10 -9 (B/ '/zG) 2 (n/cm 3 ) . This small value suggests that 
multiple strea ms play an important role in the dynamics of shock and solitary waves in the interstellar medium. 

Simulations |ll|Jl2f| have suggested that multi-stream situations are unstable, and thus would require a time- 
dependent treatment. We take an analytical approach to the question of stability by applying the collisionless 
Boltzmann equation, 

| + v .|Z + e(E+ i xB ).|^, a, 

to both species comprising a pair plasma. Since a<l, the initial fields have negligible impact on length scales of the 
order of c/u> p , where ui p is the plasma frequency. Then we could analyze the stability locally (over distances < c/ui p ) 
with the fields completely neglected. Their primary role is to serve as catalysts which induce multiple streams. With 
this, Eq. (jl]) becomes 

^+v.^=0. (2) 
at ax 

The initial two-stream distribution function satisfying Eq. (||) is locally of the form 

71 

fa[x Q , p) = - [6 (p x - p x0 ) + S(p x + p xQ )} S (p y ± p y0 ) 5 (p z ) . (3) 

The two species have values of p y o that are equal in magnitude and opposite in sign because the magnetic field bends 
particles of opposite charge in opposite directions. As a consequence, it is impossible to Lorentz boost to a frame 
where both species have vanishing p y o . It is this nonzero value of p y o that distinguishes our model from the classical 
two-stream instability. 

We now calculate the effects of a perturbation in the fields and the distribution functions with the spacetimc 
dependence exp i(k • x — uit). Faraday's Law implies that the perturbed magnetic field is related to the perturbed 
electric field by k x E = wB. Using this in conjunction with Eq. (|l|) while keeping only first-order terms in the 
perturbed quantities yields 

(-tw + ik • v) Sf = -e^ • [E + w~ 1 v x(kx E)] . (4) 

As suggested by fig] , we specialize to the case in which E is along the anisotropy axis, i, and k is along the y-direction. 
Expanding the triple vector product, we find that 
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For weak fields, the relation between E and the displacement vector D is linear, D a = e a pEp. From this, the 
definition of D in terms of the polarization vector, P = (D — E) , and the relation between the polarization and 
the current, dP/dt — J , we obtain the following relation between the current and the electric field: 
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A scalar equation is obtained by contracting with E a . Since the initial fields are negligible by assumption, the 
permittivity tensor may be decomposed into longitudinal and transverse parts. The transverse piece is given by 



(e a p - S a p) Ep = (e t - 1) E a 



And so we obtain an equation for e t : 
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The delta functions make the integrations straightforward. We integrate by parts to obtain 
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where the velocities have been normalized by the speed of light. To obtain a dispersion relation for transverse, 
electromagnetic waves, we set the permittivity tensor equal to k 2 /uj 2 . For simplicity we normalize the frequency in 

,, and the wave number by io v j c. Solving for to leads to a cubic equation 



terms of the relativistic plasma frequency, 
in lo 2 , 



10; 



= + (2*4 - 2 - 2fcV - k 



(k\ 



4 
2/0 



2k 4 v 2 y0 



±k 2 v 2 y0 



2k z v 



2„,2 



2k z v. 



2„,2 



xO v ijQ 



> UJ 



,6 4 
" v y0 



2fc 4 ^ Q 



u y0) 



(10) 



This cubic depends on two parameters, viz., the two components of the velocity. 

The exact solution to this equation is quite complicated and we will not pause to write it in its entirety. Instead, 
we mention a few special cases. First note that when p X Q = 0, the frequency is strictly real [See Eq. (JsJ) ] . Thus any 
instability must arise from having overlapping streams of the same particle species. The solution w (k) is plotted for 
two cases in Fig. [l]. As expected, when p y o — 0, the frequency is either strictly real or strictly imaginary. This is 
the classic two-stream instability. As p y o becomes comparable to p x o, we see that the plasma is still unstable; but in 
contrast to the previous case, both real and imaginary parts of the frequency are nonzero. The growth rate increases 
monotonically with k and asymptotically approaches a number of the order of the plasma frequency, while the real 
part of the frequency grows approximately linearly with k when ck 3> oj p . The instability is most powerful for large fc, 
or alternatively, over small distances. This result is self-consistent with our initial restriction of studying phenomena 
over distance scales shorter than c/uj p . 
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FIG. 1. The solid line depicts the growth rate of the standard two stream instability as a function of wave number. Each 
stream has a Lorentz factor of 100. The dashed and dotted lines show the real and imaginary parts, respectively, of the 
frequency for the modified two-stream instability discussed in §11. Here also the overall Lorentz factor is 100, but the x- and 
2/-components of the velocity are equal in magnitude. Only when the y-component is nonzero are both the real and imaginary 
parts of the frequency are nonzero. 

III. BASIC FLUID EQUATIONS 

In the rest frame of a planar shock front, cold material is flowing at a highly supersonic speed from x — — oo. 
The results derived in the previous section imply that a cold, low-cr plasma inevitably develops multiple streams 
and becomes unstable around the shock front. Simulations have shown that the instability tends to heat the plasma 
and eventually saturate JlT| , ^2[ . After saturation, we assume for simplicity that the collisionless plasma approaches 
a steady state. In this section we develop the time-independent equations necessary for the fluid description of this 
state. 

We consider the following boundary conditions at x — — oo. In the plasma's rest frame, we allow for a magnetic field 
in +z-direction. A planar shock wave travels along the x-direction with a constant Lorentz factor, 7^. We next boost 
to the rest frame of the shock and assume that the plasma has reached a quiescent state, so that all physical quantities 
are independent of time in this frame. Moreover, we take the shock front to be infinite in the y- and z-directions; 
therefore, all physical quantities can only depend on the ^-coordinate, which represents the perpendicular distance to 
the shock front. The Lorentz boost to the shock frame generates an electric field in the y-direction. Other components 
of the fields and velocities are assumed to be initially zero. 

For the sake of simplicity, we work strictly with dimensionless quantities. Electromagnetic fields are normalized by 
a reference value, Bq; spatial components of the four-velocity are normalized by a reference value jqPqc, while the 
time component are normalized by 70. Distances are normalized by jQm r c 2 / |e| Bq where m r is a reference mass to 
be determined by the particle species of interest. Finally, we define two dimensionless variables a = Bq / 8irjQm r noC 2 
and T = Pq / nQfn r c 2 . Note that these reference values are not, in general, equal to the values given by the boundary 
conditions at x — —00. 

The equations of continuity can be obtained by integrating Eq. (|l|) over momentum. The resulting equations can 
be immediately integrated yielding 



nu x — constant, 



(11) 



where u x is the x-component of the four-velocity. Faraday's law shows that the y- and z-components of the electric 
field are constant, and can therefore be evaluated from the boundary conditions. The condition, V • B = 0, ensures 
that the ir-component of the magnetic field remains zero. Poisson's equation determines how the longitudinal electric 
field develops, 
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The subscript e denotes quantities pertaining to electrons, and the subscript i denotes quantities pertaining to the 
ions (either positrons or protons). The x-component of Ampere's Law is satisfied identically; the y and z components 
are 
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We assume that each species of particles has a stress-energy tensor of the form T^ v = diag(p, P. P. P). Though 
P is not strictly the pressure in a collisionless fluid, we will refer to it as such throughout this paper for the sake of 
brevity. To these we add the stress- energy tensor of the electromagnetic field, and set the divergence of the sum equal 
to zero to obtain the Euler equations [see Eqs. ([l7|)-([T9|)]. Only three of these equations are independent; therefore, 
we will only concern ourselves with the spatial components. To complete our set of equations, we need two equations 
of state for the electron energy density and pressure, and two equations of state for the correponding ionic quantities. 
We assume that an adiabatic description is valid, with adiabatic index T between one and two @||,|9j. Then we 
close our set of equations with 



P = Po 



(15) 
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Eqs. (11), (plf) , and ( |l6| ) allow us to eliminate n, P, and p, and so express the remaining equations completely in terms 
of the four-velocities and electromagnetic fields. As described above, each species has three equations to describe its 
four- velocity, which are given by Eqs. ([l7|)-(19). Because of our boundary conditions on B v and u z , it is now clear 
that they will remain zero. 
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Our problem is to solve the remaining six, coupled, non-linear differential equations for B z (or hereafter simply 
B), E x , and the x and y components of the electron and ion four- velocities. This problem is simplified by the fact 
that three integrations can be immediately performed. These represent the conservation of the x and y components 
of momentum flux and the conservation of energy flux. The equation corresponding to the conservation of the z- 
component of momentum flux is satisfied identically. If T^ v represents the combined energy-momentum tensor of all 
species of particles and the fields, then these integrals may be obtained most readily by realizing that T Ml = constant. 
The resulting equations are given below. 
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IV. PAIR PLASMAS 



In this section we apply our model to pair plasmas. We choose the reference mass in the definition of a and T to 
be the electron mass. Because both species have the same mass, the x-components of their four-momentum will be 
equal, and the y-components will have equal magnitude and opposite sign. FromEq. (jH), we see that no longitudinal 
electric field will develop. Thus we need to determine spatial structure of the three quantities u x , u y , and B, which 
are linked by Eqs. (^0| ) and (|2l]). The simplest possible solution is one that is constant in space. This solution can 
be used to assemble a simple picture of a low-er collisionless shock. A plasma, with zero temperature at x = — oo, at 
some point encounters a shock wave generated by the modified two-stream instability of Sec. [UL We expect the shock 
width to be of order c/lo p ( <C |e| B/jmc), and so we approximate it as a discontinuity. Then Eqs. (p()|)-(p2|) may be 
used to obtain a constant solution for the post-shock plasma |isf| . In situations where er -C 1, the jump conditions 
approach the classical hydrodynamic limit, given by p0[ , as might have been expected on the basis of our equations 
of state. 

Next we examine under what conditions the constant solution is unique. Since the jump conditions of |[9[ were 
unique, the equations of a hot plasma cannot admit continuous shock wave solutions. Solitary wave solutions for a 
cold pair plasma were explored by Q. Here we extend their results to hot plasmas. From Eqs. ( p0| ) and (pl|), it is 
clear that all quantities must be bounded. From Eq. (|2l]), we can consider B 2 to be a function of u x . As u x goes to 
zero through positive values, B 2 — > — oo. As u x increases, B 2 increases monotonically, eventually becomes positive, 
reaches a maximum, and then decreases monotonically. Therefore the reality of B implies the existence of a maximum 
and minimum value of u x , and also an upper limit on the magnitude of the magnetic field. We do not consider the 
case where u x is less than zero so as to avoid instabilities. 

Eqs. ([T~T| ) and (18) show that the derivatives of the four- velocity diverge at a certain value of u x . This occurs 
because of the presence of the pressure terms, and is absent in the zero-temperature case. This divergence might 
be remedied by more exact equations of state. The presence of this divergence, however, complicates the case under 
consideration. This divergence occurs at a value of u x such that B is real; in fact, it occurs at value of u x such that 
the right hand side of Eq. (|2l]) is maximized. So we must now not only look for a bounded solution, but a bounded 
solution for which both turning points of u. x are on the same side of this singularity. 

On each side of the singularity, there is exactly one point where B — 0. Oscillatory solutions require at least two 
extrema of u x , so it is clear that we must find at least one extremum of u x where u v = 0. But in fact we must find 
at least two. For suppose that there exists an oscillatory solution on the right of the singularity. Then there must 
exist a minimum of u x such that u y = in between the singularity and the point where 5 = 0. Examination of 
the derivative of Eq. (^j) indicates that B > for a minimum of u x at such a point. Since the only place where B 
could go through zero is at a maximum of u x , we see that if B ever becomes less than zero there no longer exists the 
possibility of minimum, and so the 'solution' will inevitably become singular. Applying similar reasoning to the left 
side of the singularity, we see that at any maximum of u x , B > 0. We conclude as follows: any possible oscillatory 
solutions must have B > everywhere. As a corollary, there must exist at least two turning points of u x (and also 
B) where u y = 0. 



TABLE I. Minimum values of a required for oscillating, steady-state solutions in a collisionless pair plasma. A value of 4/3 
is assumed for the adiabatic index. 
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In addition to these restrictions, we can show that there cannot be any oscillatory solutions on the left of the 
the singularity. At an extremum of u x , both u y must vanish and the second derivative of u x must have the sign 
appropriate to whether the extremum is a maximum or a minimum. This sign is completely determined by whether 
or not B is greater or less than j/u x . Since we are only interested in B > 0, we can equivalently compare B 2 and 
7 2 /u 2 . Note that the right-hand side of this inequality is a monotonic decreasing function of u x , and the left hand 
side is given by Eq. (|2l| ) and is montonic increasing in the region of interest. Thus the two curves will intersect either 
zero or one time. We see that maxima could occur only to the left of this intersection, while minima could occur only 
on its right. Since there is no possibility of having a maximum greater than a minimum, we must conclude that there 
exist no oscillatory solutions on the left of the singularity. 

So for an oscillatory solution to exist, it must be on the right of the singularity, and there must be at least two 
points where u y — 0. We can find any such points by equating Eq. (pi]) and the square of Eq. ( |20| ) with u y set equal 
to zero. For a given u x , Eq. ( pl| ) gives the actual value of B. Eq. ( |20| ) gives only an upper limit since increasing it 2 
can only decrease B. Thus, if there is an oscillatory solution in some range of u x , the curve defined by the square 
of Eq. (20) with u y = must be greater than the curve defined by Eq. (|2l]) in the region between the extrema. The 
general solution of these equations, involving four parameters, a, 70, T, and T, with T not necessarily rational, cannot 
be expressed in terms of elementary functions. Even if we take the simpler cases of T = 4/3 or 3/2, corresponding to 
three or two dimensional relativistic velocity dispersion (at least in the collisional case), we still must numerically solve 
a polynomial of high degree in a three-dimensional parameter space. We therefore concentrate on the ultrarelativistic 
case. Of particular interest is whether, for a given temperature and large Lorentz factor, there is a range of a which 
yields an oscillating solution. We examined this question numerically for T = 0, 10~ 2 , 1, 10 2 and 10 4 . Our results, 
shown in Table |, indicate that oscillatory solutions are indeed possible. Similar to when T = 0, only magnetically 
dominated plasmas exhibit these solutions. The effect of a nonzero temperature is to increase the minimum value 
of a for which a soliton solution emerges. These results show that there are no oscillatory solutions for a <C 1 in a 
relativistic pair plasma. 



V. ELECTRON-PROTON PLASMAS 
A. Cold protons, hot electrons 

Next we investigate the steady-state properties of a plasma consisting of cold protons and hot electrons. In this 



section, we choose the reference mass in the definition of a and T to be the proton mass. As derived in Sec. Ill, 
there are only three independent quantities once Eqs. (|20|)-(|2^) are accounted for. As with pair plasmas, a trivial 
steady-state solution exists in which all quantities remain constant throughout space. We now consider perturbations 
about this solution. 

Because the electron-to-proton mass ratio is very small, quantities pertaining to the electrons vary considerably 
over a much shorter distance than quantities pertaining to the protons. In order to reduce the number of dependent 
variables and make analytical progress, we begin by investigating the spatial structure of the plasma on length scales 
of order the proton Larmor radius. We assume that the electrons undergo oscillatory motion about the spatially 
uniform solution, and then average the governing equations over several electron orbits. Because the motion of the 
electrons is bounded by assumption, we see that their averaged velocity derivatives vanish. Eqs. (|l7|) and (1^) for the 
protons go through essentially unchanged. The electrons' contribution to the right hand side of Eq. ([14 ) vanishes, 
and their contribution to Eq. ( |l2] ) is just a constant number. In Eqs. j20|)-(|22]) we average over powers of the electron 
bulk velocity. So as not to complicate the resulting equations, we take < u v x >—< u x > p . Our justification for doing 
this rests in the facts that p « 1 in almost all of the terms, and in that we are looking for solutions that stay close to 
the constant solution. 

With the above averaging, the problem is again reduced to solving for a single dependent variable given two fixed 
parameters, a and 70. Here a is defined in terms of the proton mass. Our averaging procedure evidently has the 
added advantage that there is no explicit dependence on the electrons' equation of state. To obtain a single equation, 
we use Eq. ( J22] ) to eliminate the electric field and Eq. ( pl| ) to eliminate u x . Eq. ( |20| ) can then be used to obtain a 
quadratic equation for it 2 with coefficients depending on the magnetic field: 
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To eliminate u y , we define a new independent variable dr — dx/u x . Then, (the averaged) Eq. fll4| ) illustrates that u y 
is simply proportional to the derivative of B. To obtain extremal values of B, we set u y = in Eq. (^3|). The result 
is a quartic polynomial which can be immediately deflated since B = 1 is a double root. The resulting quadratic can 
be solved to obtain the two additional roots 



B = -l±Poy/4 + 2/(T. (24) 



We now investigate what bounded solutions are possible. Eq. (|23|) can easily be solved with the quadratic formula 
for (dB/dr) 2 yielding 




1 \ „ „o A 1 



, 7; + 4, - — ^- J — \ /2fl» - (^2+ -j + # ^1 - -j + - (25) 

A sign ambiguity results; however, once the sign of the radical is chosen, it will remain unchanged provided that 
00 > l/v2- For lesser values of /3q, the radicand may go through zero; but since the zeros of the radicand are 
in general different from the extrema of B, this would violate the condition that u y be real and there would be 
no solution. Thus, once the sign of the radical it will remain fixed. Additional intuitive insight may be gained by 
considering B to be a generalized coordinate ||. Then the differential equation has the form ^B 2 + $ (B) — , i.e. 
Kinetic Energy + Potential Energy = 0. From this we can see that 'energy' is conserved; thus, there will be no 
shock wave solutions. 

We now investigate the possibility of oscillatory solutions for a -C 1. There are two possible extremal values of 



B (given by Eq. 24) which may be turning points of a physical solution. First note that if cr <SC 1, then the root 
B = 1 occurs when we choose the minus sign in Eq. ^5|. To determine the sign which would make the other two roots 
extrema, we substitute the values from Eq. 24 into Eq. |25|. For /3q ~ 1 and a <C 1, the terms outside of the radical 
of Eq. £5] are negative, indicating that the plus sign is the appropriate choice. We can set these terms to zero to see 
what choice of parameters results. Solving for a, we see that if /3q > l/\/2, it is impossible to have B given by Eq. |24| 
be roots of Eq. ^5] with the minus sign. Therefore, if a oscillating solution exists in the relativistic case, the extrema 
of B are given by Eq. [24| and we must choose the plus sign in Eq. [25[ 

However, the requirement that u x > precludes the possibility of an oscillating solution. For we have 

-.-i-^-U-g. m 

Provided /3q > l/\/2, the point B = 1 is in between the two extremal values of B. Therefore, we evaluate Eq. ^ for 
B = 1. This gives u x = — 1 + 4cr, which is certainly less than zero in many astrophysical situations where a ~ 10~ 9 
(See §11). Thus, there are no oscillating solutions, and the constant solution is unique. 



B. Hot Protons, Hot Electrons 

Finally, we consider the case where both the electrons and the protons are hot, i.e. they both satisfy the equations 
of state ( |i"5| ) and (pil). Again, we define a and T in terms of the proton mass. At x = — oo we take both species 
to be cold, but at some later stage we suppose that they have undergone a shock, in which both fluids were heated. 
We can then find the uniform post-shock solution as follows. Eqs. ( O ) and ( |l4| ) show that both species have equal 
velocities. Eq. ( P~7| ) shows that u y = for both species, and therefore Eq. ( Bp) indicates that the longitudinal electric 
field vanishes. Eqs. (^0|) and (^l|) may be used in conjuction with Eq. (|l8|) to determine J5, u x , and the total pressure. 
It is impossible from these considerations to determine how the total pressure is split between the electrons and the 
protons. Table || shows the downsteam values of these quantities for several choices of 70. We have found the jump 
conditions to be independent of a provided that a is less than approximately 0.1. 



TABLE II. The shock jump conditions for a low-a, electron-proton plasma for various shock Lorentz factors. Upstream of 
the shock, the plasma is assumed to be cold. But far downstream of the shock, the temperature has significantly increased. 
The jump conditions are essentially independent of a, provided a <<iC 1. The adiabatic index F has been taken to be 4/3. 



Lorentz Factor 


Bds/B us 




Ti + T e = Poi/n m p + P 0e /n m p 


l(f 


3.04 


3.5xl0~ a 


3.54 


10 3 


3.00 


3.5xl0~ 4 


16.7 


10 4 


3.0 


3.5xl0" 5 


77.4 
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Analytical study of non-constant solutions of these differential equations is much more complicated than in the 
preceeding cases. Here again, one must deal with coupled differential equations. One can try averaging over the 
electrons' orbits, but solving for u ytPro tons is no longer straightforward, since it is no longer easy to eliminate u XyPro tons- 
We have, however, numerically integrated these equations with er -C 1 and 70 3> 1, and have failed to find satisfactory 
solutions. 



VI. CONCLUSIONS 

We have examined the equations of a hot pair plasma, a plasma of hot electrons and cold protons, and finally a 
plasma in which both electrons and protons are hot. In none of these cases were continuous shock wave solutions 
found. However, we have seen that a cold plasma with a small, embedded magnetic field is subject to a version of the 
two-stream instability, which is a likely mechanism for the generation of collisionlcss shocks. The standard two-stream 
instability can only be saturated by nonlinear effects because of its aperiodic nature, indicating that the magnetic 
field may be significantly amplified fl8|| . But with this modified two-stream instability, both the real and imaginary 
parts of the frequency are nonzero, and so kinetic effects such as collisionless damping and resonance broadening may 
play an important role in the eventual saturation. The final level of magnetic field amplification is unclear. Since the 
growth rate of the instability is much higher than the cyclotron frequency, we expect this instability to be a dominant 
influence in the formation of a collisionless shock. 

We have also investigated the possible existence of soliton solutions in a hot, pair plasma. As with cold plasmas, 
we have found that such solutions exist only for a 3> 1. Thus, the only time- independent solution is one that is 
independent of the spatial coordinate. This simple structure has been observed in numerical simulations Jl5| , which 
do suggest spatially uniform solutions, with physical quantities equal to the values predicted by the jump conditions. 

We have also considered a model consisting of hot electrons and cold protons. Since we have averaged over the 
electrons' orbits, this model is essentially independent of the electrons' equations of state. In this case again we have 
found that the unique solution for a relativistic, low-cr plasma is the spatially uniform one. Finally, we obtained the 
jump conditions for the state of a hot plasma of electrons and protons that was initially cold and had undergone 
a collisionless shock. As in the other cases, we have found that almost all of the plasma energy is converted into 
pressure. The magnetic field gets amplified by a factor of ~ 3 in the shock frame, and the post-shock value of a 
is only amplified by a factor of ~ 7^. In many astrophysical situation, this leaves the fraction of energy contained 
in the magnetic field still far below the equipartition limit. Numerical simulations of relativistic collisionlcss shocks 
would be very useful in understanding whether these jump conditions are appropriate, and whether or not the adopted 
model of a quiescent, low-cr plasma is adequate. In particular, such simulations would also provide information on 
the relative thermal Lorentz factors for the electrons and protons, which is of fundamental importance for models of 
the relativistic blast wave in gamma-ray burst afterglows JL8[ . 
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